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Abstract 

We compute total cross sections for various slepton pair production reactions 
^Li'L-, ^L^L, ^R^R and DlDl in next-to-leading order QCD. For pp collisions at 
"v/S* = 2 TeV, we find leading order cross sections to be enhanced by typically 
35% to 40%. For pp collisions at = 14 TeV, the enhancement ranges from 
25% to 35% depending on the mass of the sleptons. We comment upon the 
phenomenological implications of these results. 



PACS numbers: 14.80.Ly, 12.38.Bx, 13.85. Qk 
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I. INTRODUCTION 



Particle physics models that incorporate weak scale supersymmetry (SUSY) are certainly 
among the most compelling possibilities for physics beyond the Standard Model (SM) [|l|. 
Numerous searches for the sparticles of supersymmetric models have been made at both 
lepton and hadron colliders. The lack of any convincing supersymmetric signal beyond 
Standard Model backgrounds has led to a variety of limits on sparticle masses . 

The best search limits on strongly interacting superpartners, the squarks q and gluinos g, 
have been obtained from the CDF and D0 experiments at the Fermilab Tevatron pp collider. 
By searching for multijet+^y events, limits of nig ^ 144 — 154 GeV (for niq » nig) and 
niq > 212 — 224 GeV (for niq ~ nig) have been found Recently, a next-to-leading 

order computation of gg, gq and qq production has been completed [Q. As a result of this 
calculation, the theoretical uncertainty due to choice of renormalization and factorization 
scales has been reduced. In addition, for a scale choice of the order of the average mass 
of the produced sparticles, the pair production cross sections were found to increase by a 
factor of 1 to 2 depending upon the subprocesses and collider. This has led to 10-30 GeV 
improvements in the mass limits coming from the CDF and D0 collaborations. 

For weakly interacting sparticles such as the charginos Wi and sleptons £r, £l and ugL 
{£ = e, /i or r), the best search limits come from LEP2 experiments. Currently, "Jti^^ > 85.5 
GeV for a gaugino-like chargino [Q. For charged sleptons, mg^ ^ 75 GeV, m^^ > 59 GeV 
and nifj^ > 53 GeV (where the exact bounds depend on the lightest neutralino mass ni-^^) 
1^. In addition, sneutrinos are required to have mass nic, ^ 41.8 GeV from measurements 
of the invisible Z width |^ (three degenerate generations of invisibly decaying sneutrinos 
have been assumed). Search limits from the Tevatron collider on charginos and neutralinos 
can be obtained by looking for clean trilepton events from W1Z2 —^3i + Ifj, production. 
Expected event rates are very model dependent, so that limits are best placed as contours 
in parameter space of specific models 0. 

The possibility of searching for slepton production at hadron colliders has been studied 
in Refs |1T0|JTT[] . In gravity- mediated models of supersymmetry breaking (SUGRA models) at 
Tevatron conditions, direct slepton pair production occurs at low rates, although it may be 
possible to extract a signal at one of the upgrade options such as the Main Injector, which 
will operate at = 2 TeV and accumulate ~ 2 fb~^ of integrated luminosity, or TeV33, 
which will operate at the same energy, but amass ~ 25 fb~^ of integrated luminosity |]1T 



For models where SUSY breaking is communicated by gauge interactions |jT2|, slepton pair 
production can lead to final states with two isolated photons plus two isolated leptons, a 
signature for which SM backgrounds are expected to be small |1^ 



For the CERN LHC operating at = 14 TeV with projected data sets of 10 - 100 
fb~^ of integrated luminosity, SUGRA models predict the possibility of detection of directly 
produced slepton pairs [lTll , p!4| . This is likely possible in gauge-mediated models as well. 
SUGRA models predict slepton masses of 100 — 400 GeV ought to be detectable; this mass 
range also corresponds to the most favored range of slepton masses expected if the lightest 
neutralino makes up the bulk of dark matter in the universe ||15|| . 

Direct slepton pair event rates will be low at an upgraded Tevatron, while at the LHC, 
slepton reconstruction will be difficult. Theoretical improvements of the Born level analy- 
ses |lTO|JTl[] of slepton pair production are essential. In this paper, we present the leading 
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Quantum Chromodynamics (QCD) corrections to the direct production of slepton pairs at 
hadron colliders. The calculation is described and formulae are given in Sec. II. Numerical 
results are shown for the Fermilab Tevatron and CERN LHC collider in Sec. III. The main 
result is that the slepton pair production cross sections are enhanced by 35-40% for the 
Fermilab Tevatron collider, and 25-35% for the CERN LHC collider relative to the leading 
order results. The theoretical uncertainties are addressed and it is found that the total cross 
section results should be accurate to ±15%. Finally, we comment upon the implications of 
these results for slepton search experiments. 



II. ANALYTICAL RESULTS 



In this section we briefly describe the calculation of QCD corrections to the direct pro- 
duction of slepton pairs at hadron colliders and then give the corresponding cross section 
formulae. The results are very similar to those of the leading QCD corrections []TB[ to the 



SM Drell-Yan [|T^ process. The lowest order cross sections have been discussed in Refs. 
and 



The colorless nature of the slepton pair final state allows one to factorize the production 
cross section into hadronic and SUSY pieces. The SUSY piece corresponds to the decay of 
a virtual SM gauge boson of mass into a pair of slepton particles with masses mi and 
m2- The decays are 



Z* (or 7*) 
Z* (or 7*) 
Z* 



Ul , 
IrIr, and 



(2.1) 



Using the gauge invariance of the hadronic piece and the slepton Feynman rules, from Ref. 
p!9| for example, the two body final state phase space integrals in the SUSY piece may be 
performed analytically. The result is proportional to (3^ where 



/3 



2(m\ + m\) 



+ 



1/2 



(2.2) 



Having calculated the SUSY piece, the QCD corrections to the hadronic piece are stan- 
dard and correspond to the production of a virtual SM gauge boson. Gluon-bremsstrahlung 
is added to the leading order processes, giving the subprocesses 



qq W* g, and 
qq Z* (or7*)5(. 



(2.3) 



Singularities corresponding to soft gluon emission cancel when the ultra-violet renormal- 
ized virtual diagram contributions are added. The remaining singularities are initial state 
collinear in origin and are mass factorized into scale dependent parton distribution functions. 
We use the modified minimal subtraction (MS) scheme throughout. 
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At this order, in addition to the annihilation processes ( |2.3| ) there are also contributions 
from Compton scattering diagrams that may be obtained from those of ( p.3|) by crossing. 
There are no soft singularities for these processes, only initial state coUinear singularities. 

The total production cross section is given by a convolution of scale dependent parton 
distribution functions and hard scattering partonic subprocess cross sections. 



(^=[ ^^'^^ L ^2 C^Xsi ^ fi/A{xA,f^f)fj/B{xB,f^f)^^ (2.4) 

J(mi+m2)2 [ij=q,q 



+ H [fi/A{xA,f^f)fg/BixB,f^f) + fg/A{xA,f^f)fi/B{xB,fif) 
i=q,q 



where fi/A{x,fif) are the parton distribution functions for parton i in hadron A with mo- 
mentum fraction x evaluated at the factorization scale /i/. S* is the hadron-hadron center of 
mass energy squared which is related to s, the parton-parton center of mass energy squared, 
via s = xaXbS. Defining z = Q'^/s the hard scattering partonic subprocess cross sections 
are given by 



da, 



qq 



dQ^' 




(2.5) 



and 



da, 



qg 



dQ^ 



ao- 



as(/ir) 1 

27T 2 



+ z--z' + 2Pqq{z)[\n 



'l-z) 



(2.6) 



where cxo are given below, fir is the renormalization scale, and Pqq and Pqg are the splitting 
kernels PO] given by 



Pqqi^) 



Pqg{^) 



1+^2 3 ^ ,^ 



z' + il-zf 



(2.7) 



For the W exchange process. 



5767r 2 

with 



^, = ^^z\Dw{Q^)\' (2.8) 
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DxiQ') = , ,^ . (2.9) 



1 

-Mf+tTxMx 

For the 'j — Z exchange processes, 



W 



with 



J'.ziQ') = \Dz{Q')\\L'^ + Rt)4 + J ^ (2.11) 



+2(1, + fi,)c,i£^!!^|£!£*£i«!i^|c,(Q^)p 

9 ^ 



The constants in J^^z are given by 



Ti = ±1/2 (weak isospin), (2.12) 
Cf = 2/3, —1/3, —1, (electric charge), 
Lf = 2Ti - 2e/ sin^ 9w, 
Rf = — 2e/ sin^ 9]y, 



and depending on whether one has a left- or right- SUSY partner, 

= Lf or Rf , 

for slepton /. 

III. NUMERICAL RESULTS AND IMPLICATIONS FOR DETECTION 

Using the results described in the previous section a program was constructed to calculate 
slepton pair production cross sections according to Eq. (p.4|) . For input, standard values 
of the electroweak parameters were used along with the CTEQ4M parton distribution 
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functions ||2T|. A two-loop expression for as was used with fixed nj = 5. The value of ^qcd 
was taken from the parton distribution set. 

The next-to-leading order (NLO) cross section results are presented in Fig. 1 for a) pp 
collisions at = 2 TeV and h) pp collisions at = 14 TeV as a function of slepton 
mass. For ^Li'L production, we have assumed m^^ = m,>^. Then for both the Tevatron and 

the LHC, we see that W* ih^L is the dominant production mechanism, while £_r£r is the 
smallest. The Il^l and z/^z^l cross sections are intermediate between these two cases, and 
of comparable magnitude to each other. 

The results of Fig. 1 were produced by setting the renormalization scale /ir equal to the 
factorization scale /i^; the square of the scale was taken to be the invariant mass of the 
slepton pair . To get some idea of the size of the uncalculated higher order terms we 
varied the scale from Q/2 to 2Q. The results changed by ±5% for the Tevatron and ±2% 
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for the LHC The largest uncertainty comes from the parton distribution functions (PDFs). 
For a crude estimate we use a set that assumes a different functional form for the gluon 



and sea quark at the initial evolution scale, the MRS Set Rl |]22|. This leads to an increase 



in cross section of approximately 10% for both Tevatron and LHC. We therefore cautiously 
estimate the theoretical uncertainty to be ±15%. 

In Fig. 2, we plot the ratio of next-to-leading order cross sections to leading order cross 
sections, again versus slepton mass for the Tevatron and LHC colliders. The leading order 
results were computed using CTEQ4L parton distribution functions. 

For the Tevatron collider, the various cross sections have an enhancement of 36 — 39% 
depending upon the slepton mass. The correction is nearly the same for all four slepton 
pair production reactions, as might be expected by the form of the correction given in Eqns. 
Q and (^). 

The most promising detection strategy for direct production of sleptons at the Tevatron 
in SUGRA-type models is via detection of acollinear dilepton pairs from, for instance, pp 
£r£r ii + Ifrp. Detailed simulations of this reaction have been performed in Ref. ||11|| . 
The major SM background to the clean dilepton signal comes from WW — > ^^' + I^j.. 
Ref. evaluated two case studies for slepton pair production. After cuts, case 1, with 
~ 80 GeV, had 18 fb of signal, while case 2 with m^^ ~ 102 GeV had 9 fb of signal. The 
background from WW at lowest order was evaluated at 36 fb. Thus, with 2 fb^^ of integrated 
luminosity, neither signal case would be visible by itself at the So" level. Incorporation of 
the NLO enhancement to both signal and background {K ~ 1.3 for WW production at the 



Tevatron ^3[) pushes case 1 above the 5cr level, and into observability. In general, for the 
Tevatron, supersymmetric clean dilepton pair events can also come from chargino/neutralino 
pair production, and the net signal would likely come from a variety of different sparticle 
production reactions. The reach for SUGRA type SUSY models in the clean dilepton channel 
has been presented in Ref. 

For the case of the CERN LHC collider. Fig. 2h shows that the various slepton pair 
cross sections have an enhancement of 25 — 36% depending upon slepton mass. In Ref. 



TT| , P^ , it was found that a clean ^^ + I^j, signal could be seen above SM backgrounds for 
100 ^ < 400 GeV. With hard enough cuts, the signal was shown to come almost 
exclusively from slepton pair production rather than from other SUSY particle production 
processes. However, direct reconstruction of parent slepton masses would be extremely 
difficult. Possibly the best option for a slepton mass measurement at the LHC might be 
from matching the total signal rate to expected rates from simulations of slepton production 
and decay for different mass sleptons. In this case, the most precise estimate of the slepton 
production cross sections would be needed. 



IV. CONCLUSIONS 

The total cross sections for various slepton pair production reactions ^li^l, ^l^l, ^h^/j and 
i^L^L were computed in next-to-leading order QCD. For pp collisions at \fS = 2 TeV, leading 
order cross sections were enhanced by typically 35% to 40%. For pp collisions at y/S = 14 
TeV, the enhancement ranges from 25% to 35% depending on the mass of the sleptons. The 
theoretical uncertainty resulting from variations in the scale and parton distributions was 
found to be approximately ±15%. The NLO enhancements of the cross sections at Tevatron 
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energies push some predictions for signals at the MI above the ba leveL At the LHC. where 
a slcpton mass measurement may be made indirectly on the basis of event rates rather than 
by reconstruction, the NLO predictions for slepton pair production rates would be essential. 
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Figure Captions 



Fig. 1 The NLO total cross section for pair production of a single generation of sleptons: 
cli^l + Cit'L, clCl, eRtR and vl^Li via the Drell-Yan mechanism versus slepton mass, 
for a) pp collisions at = 2 TeV and b) pp collisions at VS — 14 TeV. The slepton 
masses are assumed to be degenerate. We convolute with CTEQ4M parton distribution 
functions. 

Fig. 2 The ratio of NLO to LO slepton pair production cross sections versus slepton mass, 
for a) pp collisions at = 2 TcV and b) pp collisions at = 14 TeV. The slepton 
masses are assumed to be degenerate. Wc convolute with CTEQ4M PDFs in the NLO 
case and CTEQ4L PDFs in the LO case. 
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